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(54) Multiuser detection in CDMA using a correlation matrix 



(57) A radio communications receiver which oper- 
ates to detect radio signals and to recover data repre- 
sentative of the radio signals in the presence of contem- 
poraneously detected interfering signals, the radio com- 
munications receiver comprising means to detect the ra- 
dio signals (6, 8) and to generate digital base band sig- 
nals representative of the radio signals, means to form 
a correlation matrix (12) from the base band signals, a 
data processor (16) which operates to generate a factor 



(R) of the correlation matrix (S), and to recover the data 
from the factor, wherein the data processor (16) oper- 
ates to generate the factor (R) by calculating first ele- 
ments of the factor (R) from the correlation matrix (S) 
using a factorisation algorithm, in accordance with a pre- 
determined structure of the factor (R), and forming sec- 
ond elements of the matrix (R) by selectively copying 
first elements in accordance with the predetermined 
structure. 
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Description 

[0001] The present invention relates to radio communications receivers and in particular to radio communications, 
receivers which operate to detect radio signals and to recover data representative of the radio signals in the presence 
5 of interfering signals. The present invention also relates to methods of receiving radio signals and recovering data 
representative of the radio signals in the presence of interfering signals. 

[0002] In data and signal processing applications, there is often a requirement to factorise a matrix of signal samples 
to the effect of providing a factor matrix having particular characteristics. For example, the desired matrix factor may 
be upper triangular or lower triangular in that the lower or upper part of the matrix may be filled with zeros. As a result, 
10 the factorised matrix may be interpreted as a set of simultaneous equations and these equations solved in a manner 
facilitated by the upper or lower triangular form. 

[0003] An example where such matrix factorisation is required is in radio communications receivers, and in particular 
radio communications receivers which are arranged to detect radio signals and generate data from these radio signals 
in the presence of contemporaneously detected interfering signals. Such is the case with radio systems arranged to 

is operate with code division multiple access. Code division multiple access is a radio communications technique in which 
data to be communicated is combined with a spreading code in some way before being modulated onto a radio fre- 
quency carrier signal. At the receiver the radio signals are detected and the data is recovered by de-spreading the 
radio signals with knowledge of the spreading code used at the transmitter to form the radio signals. As a result of the 
spreading of the data, the receiver is able to generate a gain in the power of the detected signal with respect to noise 

20 and other signals for the wanted signal. As such, signals from a plurality of transmitters may be contemporaneously 
transmitted using different spreading codes and separated at corresponding receivers to the effect that the data which 
the radio signals represent may be recovered in spite of the presence of the interfering signals from the other trans- 
mitters. However in order to detect the data from a wanted radio signal, the receiver is required to reject the unwanted 
signals which are contemporaneously detected with the wanted signal. As such, the receivers require extensive signal 

25 processing capabilities in order to effect the task of detecting the data in the presence of the unwanted signals. This 
represents a task of some considerable complexity. As such, it may not be possible to detect the data in real time, 
because state of the art signal processors are not able to execute the number of calculations required for the signal 
processing algorithms before the data must be presented at an output. However, it is not always necessary to solve 
signal processing problems exactly. It is known to be possible to only compute an approximation to the exact solution 

30 without degrading the overall performance of the system. Computing an approximation to the solution reduces the 
amount of time spent solving the problem, reduces the amount of hardware required for implementing the algorithm, 
reduces the power consumption of the device and in some cases makes it possible to achieve real-time behaviour in 
the first place. 

[0004] A part of a process most appropriate for detecting data communicated in accordance with a code division 
35 multiple access system requires the factorisation or decomposition of a correlation matrix. In order to effect real-time 
operation with known signal processors, it is necessary to provide an approximation to an exact solution of a factori- 
sation of this correlation matrix. In an article entitled "Real-time feasibility of joint detection CDMA" by J. Mayer, J. 
Schlee and T. Weber, published in the proceedings of the second European Personal Mobile Communications Con- 
ference, Bonn, Germany, September 1997, pages 245 to 252, an approximation of a matrix factorisation is described 
4> using a known numerical approximation process. However, the computational effort needed to perform this multi-user 
detection is still considerable. 

[0005] A technical problem therefore exists in further reducing the computational complexity, or correspondingly 
increasing the accuracy which can be achieved using an approximation method to generate factor of a matrix. 
[0006] The technical problem is addressed generally by considering the generation of the matrix factor by analysing 
45 the structure matrix factor and determining from the structure which of the elements of the matrix must be calculated 
ex actly u sing a factorisation algorithm, and which of the elements may be copied from the components calculated 

[0007] According to the present invention there is provided a receiver which operate to detect radio signals and to 
recoverydata representative of the radio signals in the presence of interfering signals, said radio communications re- 
[ 50 ceiver/comprising means to detect said radio signals and to generate digital base-band signals representative of the 
radip signals, means to form a correlation matrix from the base-band signals, and a data processor which operates to 
emerate a factor of the correlation matrix and to recover the data from the factor, wherein the data processor operates 
"to generate the factor by calculating first elements of the factor from the correlation matrix using a factorisation algorith m 
in accordance with a pre-determined structure of the fapipijnatniand forming second elements of th e matrix by se- 
55 lectively copying said first elements from a first row irfaccordance with the pre-determined structure. 

[0008] The pre-determined structure of the factor matrix may be determined in accordance with a band structure 

i form of the calculation and the technique used to.represent.the^xlata-as radio signals. 
^[0009] \ The first and second elements maybe rows or columns oUKern^Wx^The term rows as used herein may be 
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substituted for the term column in that by invertin^either the correlation matrix or the factojjnatrix 6r bot)v/the columns 
become rows and the rows become columns. Thprfactor may be and uppenriangular matrix or a lower.tri^ngular matrix. 
The radio signals and the interfering signMs-rriay be generated in accordance with code division multiple access. The 
code division multiple access may be time division-code division multiple access. 
5 [0010] The known approximation methods for factoring a correlation matrix as disclosed in the above mentioned 
reference entitled, "Real-time feasibility of joint detecti on CDMA", by Mayer, Schtee, Weber, have been guided by a 
mathematical view of the calculation and the treatment has been restricted to one specific algorithm for computing this 
part, namely the Cholesky algorithm. The present invention, however results from an algorithmic analysis of the matrix 
factorisation, from which a significant reduction in complexity is facilitated. 
10 [0011] For the particular example of signal processing in a CDMA receiver, an approximation method greatly increas- 
es the capability of the receive r to cope with a technical problem known to those skilled in the art as the near/far 
problem, where different users are received with very different signal powers. As contemporaneously received unwant- 
ed signals with comparatively high signal powers tend to prevent the detection of wanted radio signals having a com- 
paratively low signal power. 

15 [0012] According to an aspect of the present invention there is provided a method of detecting a radio s ignal and 
recovering dat a representative of said radio signal in the presence of interfering signals, said method comprising the 
steps of; 

- detecting said radio signals and generating d igital base band sig nals repr esentative of said radio signals, 
20 . generating a correlation matrix (S) fr om said baseband signals, 

- generating a matrix factor (R) of said correlation matrix, and 

recovering said data from said factor (R), wherein said factor (R) is generated in accordance with the following 
operations; 

calculating first elements of the factor (R) from the correlation matrix (S) using a factorisation algorithm, in accord- 
25 ance with a pre-determiried structure of the factor (R) ; and 

forming second elements of the matrix (R) by selectively copying elements from said first rows in accordance with 
said pre-determined structure. 

[0013] According to a further aspect of the present invention, there is provided a method of generating at least one 
3a factor (R) of a matrix (S), said factor (R) being upper or lower triangular in form, said method comprising the steps of; 

calculating first elements of the factor (R) from the matrix (S) using a factorisation algorithm, in accordance with 
a pre-determined structure; and 

- forming second elements of the matrix (R) by s electively copyin g elements from said first elements in accordance 
35 with said pre-determined structure. 

[0014] One embodiment of the present invention will now be described by way of example only with reference to the 
accompanying drawings wherein; 

FIGURE 1 is a schematic block diagram of a mobile radio telecommunications system; 

FIGURE 2 is a schematic block diagram of a receiver which operates within the mobile radio communication system 
shown in Figure 1 ; 

45 FIGURE 3 is an illustrative representation of a known method of generating an approximation of a factor matrix R, 
from a correlation matrix S; 



50 



FIGURE 4 is a representation of a process of generating a factor matrix R. from a correlation matrix S according to 
the present invention; and 

FIGURE 5 is a further illustration of a process generating a factor matrix of a correlation matrix S according to the 
present invention. 



[0015] One embodiment of the present invention will now be described with reference to a mobile radio communi- 
55 cations system which operates in accordance with code division multiple access. The present invention may operate 
with time division-code division multiple access (TD-CDMA) or wide band code division multiple access (W-COMA). 
Both multiple access techniques will benefit from exploiting opportunities for approximation. The TD-CDMA system is 
required to perform joint detection of all users that are active in one time/frequency slot. Compared to a TD-CDMA 
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system, the number of simultaneously active users in a W-CDMA system is much higher. Therefore, it is currently not 
possible to jointly detect all users. However, it is possible to perform multi-user detection, where only a subset of all 
users is included in the joint detection problem and the rest is treated as interference. Approximation methods increase 
the number of users that can be included in one detection step and can thus increase the resistance of the system to 

5 the near/far effect. This in turn leads to reduced demands on the power control. 

[0016] The illustrative embodiment if the present invention will be described with reference to a TD-CDMA system. 
An example of part of a mobile radio telephone system which operates in accordance with a TD-CDMA system is 
illustrated in Figure 1. In Figure 1, three base stations BS are shown to be interconnected via a mobile network infra- 
structure, NET. Data is communicated between mobile stations MS and the base stations BS, by transmitting and 

io receiving radio signals 1 , between the base and mobile stations operating within a radio coverage area effected by 
each of the base stations BS. The radio coverage area is shown illustrated as a broken line 2, and serves to indicate 
a boundary within which radio communications can be effected with the mobile stations MS. In the present illustrative 
embodiment the mobile stations MS communicate with the base stations 8S in accordance with a TD-CDMA system. 
A more detailed explanation of how data is communicated using a time division CDMA system is provided in an article, 

15 entitled ..Performance of a Cellular Hybrid C-TDMA Mobile Radio System Applying Joint Detection and Coherent Re- 
ceiver Antenna Diversity" by 6. Blanz, A. Klein. M. NaBhan and A. Steil published in the IEEE Journal on Selected 
Areas in Communications, Volume 12, no. 4, May 1994 at page 568, the content of which is incorporated herein by 
reference. However for the present explanation it should be noted that communications using TD-CDMA systems is 
characterised in that radio frequency carrier signals via which communication of data is effected are divided into a 

20 plurality of time-slots. Each of these time-slots is assigned to a plurality of mobile stations which operate to communicate 
radio signals in the time-slot. In order to separate data communicated by the plurality of mobiles assigned to the same 
time-slot, each mobile is provided with a user specific code which is convolved with the data to be communicated. A 
receiver of the radio signals operates to separate data communicated by other mobile stations MS contemporaneously, 
in accordance with the user specific code. As will be appreciated, a simitar arrangement is effected for data commu- 

25 nicated from the base stations to the mobile stations. 

[0017] A receiver which operates within the mobile radio communications system shown in Figure 1 , is illustrated in 
Figure 2. In Figure 2 a wanted radio signal transmitted from a mobile station is represented by the line 1 , whereas 
unwanted signals transmitted by other mobile stations are represented by the lines 4. The wanted and unwanted radio 
signals are detected by an antenna 6 and fed to a radio frequency down converter mean s 8^of-the.r eceiver 1 0. The 

30 radio frequency down converter means 6 operates to c onvert the received radio signals into a base band /epresentation 
w hich are the reafter_analogue to digital converted and fed to a first data processor 12. The first processor 12 operates 
generate a correlatio n ma trix S^ of the received ra dio signals in combination with a plurality of CDMA-s preading 
codes fed from a data store 14. T he plurality of spreading codes held in data store 14 correspond to each of the 
spreading codes used to spread the data modulated onto the wanted and unwa nted radio signals. The data store 14 

35 may also include a set of training sequences w hich are embedded in the unwanted and wanted signals from which a 
channel impulse response of the communications channel through which the wanted a nd unwanted radio signals have 
passed may be estimated. The correlation matrix S t is thereafter passed to a further data processor 16 which operates 
to generate a f actorised matrix R from the correlation matrix S from which data representative d of the wanted radio 
signals is generated an output on a conductor 18. Operation of the data processor will now be described in more detail. 



40 [001 8] The present example embodiment will be illustrated with regard to a receiver 1 0, which operates in accordance 
with TD-CDMA. As with W-CDMA, the problem can be reduced to solving a system of linear equations, after the 
correlation matrix S has been generated by the receiver A rec eived radio signal which has been generated in accord- 
ance with a TD-CDMA system can be modelled with an over-determined linear equation system, as illustrated'in equfe- 



[0019] The system matrix A contains the combined channel impulse responses and spreading codes of each of the 
so signals contemporaneously received in each time slot. The vector e contains the received antenna samples and the 
vector d denotes the unknown data symbols. This system of equations is solved according to a least squares criterion. 
This may be effected, by factoring a correlation matrix, formed from the matrix A, by multiplying the matrix with a 
Hermittian transpose of itself. This is illustrated by equation (2): 




tion (1): 



45 



Ad = e 



(1) 
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A H Ad = A H e 
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[0020] The data is recovered from equation (2), by forming a factor matrix R, as determined by equation (3): 

A H A=tf*R (3) 
[0021] Where in equation (3), R is upper triangular. From equations (2) and (3), equation (4) may be formed: 

R M Rd = A H e (4) 

[0022] Equation (4) can be trivially solved by two substitution operations to recover the data vector d. To facilitate 
> the-following explanation, the correlation matrix to be factored is designated S, where S = R. 
[0023]V A key step in this process of recovering the data ve ctor d , is generating th e factor matrix R. This can be 
generated exatflyjjsjng an.algorithnxknown to those skilled in the art as the Cholesky algorithm, a description of which 
may be/found in a publication entitled "Coding, Modulation and Multiuser Decoding for DS-CDMA Systems", PhD thesis 
_byJ>Ottosson, Chalmers University of Technology, GOteborg, Sweden, 1997. Having reg ard to the number of elemen ts 
present in the/fnatrix s ))which is det ermined l y the numbe r of data symbols in the vector d, in combination with the 
length of the s preading sequen ce and the d elay spread of the channej ^pjjjse^response, the elements of the factor 
matrixT? c afrno t be calcula ted.exactly using the Cholesky algorithm, as this would either prevent real-time operation, 
or require a data processor with a processing power not cu rrent available within the state of the art. However the matrix 
S is known to hav e a structure which is Hermitian, posi tive definite an d is als o a block Toeplitz matrix with a band 



structure. As such a known method for ge neratin g an a p proximati on to the factor matrix R, is disclosed in the above 
refere nced publication entit led "Realtime feasibility of joint detection CDMA" by Mayer, Schlee and Weber. This ap- 
proximation will now be described with reference to Figure 3. 

[0024] Figure 3 provides an illustrative representation of an approximation to the matrices S and R. For the correlation 
matrix S, and the matrix factor R t the following conventions are adopted: 

S = (s fj ) t S=A H A 

R = (rtf, with R upper-triangular and R H R = S, and for both S and R t 1 zizn , 1 zjzn . 

[0025] Let the blocks of S and R be of size kxk. Then the block-Toeplitz structure of S and the approximate block- 
Toeplitz structure of R can be described as expressed in equation (5): 

S U s S «('M and r tj a r UQ+ ( 5 > 

[0026] Let S be a square sub-matrix from the upper left comer of S, which is shown in Figure 3, bounded by a dashed 
line 20. These elements of the matrix S, are bounded by an index ft, as shown in Figure 3, and described mathematically 
by equation (6): ' 1 * * "" 

S \£i£n <n , l£ j sn <n (6) 



[0027] The Cholesky factor X of S is then a sub-matrix from the upper left comer of R bounded by a dashed line 22, 
as expressed by equation (7): 



R = (ry), \Zi*n<n, \£jzH<n (7) 



[0028] The exact matrix R can thus be approximated by computing the Cholesky factor from a smaller sub-matrix 
$ of S and then filling the rest of R with copies of ele ments of R , as indicated by (5). As shown with the dashed line 
20, 22, in Figure 3, $ and l^afeplaced at the upper left comer of S and R, respectively. After computing R , shown 



5 



35 



EP 0 971 485 A1 

as darker blocks on the right hand side of Figure 3, the remaining. blocks of R, which are shown in a lighter shade, are 
copied from corresponding columns of the sub-matrix , in accordance with a pre-determined mathematical structure. 
[0029] As described in the above referencac^pbblication "Real-time feasibility of joint detection CDMA", by Mayer, 
Schlee and Weber, the actual computation of jLx\ be done with the Chotesky algorithm, which requires one square- 
5 root operation and one division for each corapdted row of $ . Let v be the degree of the block-band structure of R, that 
is, that there are only vkxk non-zero blocks directly on the diagonal and to the right of R in each block-row. To capture 
all the information of R, one needs to choose an X that has at least a size of v x v blocks, as indicated in Figure 3. 

This results in kv square root operations and kv divisions as required by the Cholesky algorithm. ^ ^\ 

[0030] This known technique for reducing a number of calculations required to generate the factor matrix Rf is base ] 
10 on a mathematical assessment of the process for generating the'matrix factor. However, the Applicant-has observed ' 
that by analysing further the structure of the factor-rqatrix Rj it is possible' toachieve further reductions in a number of 
calculations required to produce the factor matrix ftJIh-this way, the approximate block-Toeplitz structure of R is ex- 
ploited, to the effect that only a sm allpart offc Is-compute d, and th en the remaining parts filled with suitable copies of 
the co mputed pa rt. Exploiting this opportunity to reduce th e number of calculations is an im portant part of being able 
*5 to perform the j oint-detection in real-tim e. The algorithms that can be used to com pute H can be written in a form so 
that the elements of # are produced row-wise from top to botto m and within each row from left to right. While executing 
the algorithm, it is possible to leave duUhe computation of more of the elements which are substituted with copies 
from those elements calculated exactly using the factoring algorithm, according to a pre-determined approximation for 
the structure of R . This approach leads to a close intertwining of the actual execution of the algorithm and the approx- 
20 ^fnaQon process and thus results in great flexibility for controlling the approximation process. 

[00311 This general principle can be applied to all algorithms that compute elements of a factor matrix R , using known 
factoring algorithms such as the Cholesky algorithm or other known algorithms such as the Block-Schur algorithm, 
which is known to those skilled in the art Furthermore, it is beneficial for a wide range of signal processing tasks, tike 
jointrdetection in TO-CDMA systems and multi-user detection in W-CDMA systems. 

032] The pre-determined structure, from which those elements of R , which need to be calculated using the factoring 
algorithm, and those elements which may be copied, can be substantially arbitrary and can vary dynamically in ac- 
cordance with, for example, the prevailing conditions for data communications using a CDMA radio communications 
system . 

[0033] There are many forms which the pre-determined structure may take. By carefully choosing one, it is possible 
30 to only compute the most "important" elements. Copying an element according to (5) makes only really sense when 
the copied-from element exists. This means that at least the f irst k rows of R should be comput ed, because there are 
no elements from which they could be copied from. As a result of the pre-determined stru cture, which may, as in 
accordance with the present embodiment, identify a band-structure of the correlation matrix 5, only the first vk elements 
of each of k rows are non-zero. In mathematical terms, this can be expressed as shown in equation (8): 



shall we compute r^7 s js / £ * and jZvk? (8) 



[0034] Figure 4 provides an illustration of how the p re-determined structure of ft is used to determine which elements 

40 can be copied and which must be calculated using the f actoring alg orithm. As illustrated, the darker shaded blocks 24, 
off ?, on the right hand side of F igure 4 are computed and the lighter blocks 26, are copied from them. Starting from 
this rule, the region of the computed elements can easily be enlarged. This generalisation has been achieved by 
examining the process of computing the Cholesky factor from an algorithmic point of view. The state of the art teaches 
the skilled person on l y a mathematica l, high level analysis of the calculation of the matrix factor R. This mathematical 

45 view only provided knowledge of matrix-level computations. Analysing an implementation of the calculation, provides 
an insight into the structure of t he factor matrix , from which a reduction in the number of c omputations can be effect ed. 
As will be clear to those skilled in the art, the illustrative embodiment of the present invention which finds an approxi- 
mation to the generation of a Cholesky factor/?, is one example of a general approach in which analysis of the structure 
of the factor m atrix reveals a way in which a reductio n in a number of computations can be achieved in calculating an 

50 approximation to a matrix factor. 

[0035] Simulations have shown that exploiting these new opportunities for saving computational effort lead to an 
improvement over the existing solutions. These simulations have been carried out with a TD-CDMA model and have 
shown that the approximation to R according to the embodiment of the present Invention ] produces rib reduction in the 
bit error rate of the recovered jata, w ith a reduced amount of calculation s. Only the first v blocks of t he-first block-row 

55 need to be computed, th e restoT gcSh be copiedjromjhem. Com pared to the known.metbod of com puting R. the 
'embodiment reduces the amount of square root operations and divisions from vk to k. For typical channel models, this 
represents a reduction in the amount of calculations by around one third. 



6 



i 




EP0 971 485 A1 




[0036] The new approach is applicable to a wide range of algorithms, such as the original Cholesky algorithm or the 
improved Btock-Schur algorithm. It can be used in many signal processing tasks that require a matrix decomposition 
to reduce time and space complexity. For example, it can be used in the joint-detection process in a TO-CDMA system 
and for multi-user detection in W-CDMA systems. It not only provides a finer control over the amount and specific 

5 placement of elements which are calculated exactly, it also provides a means for interleaving computation and approx- 
imation. This enables the matrix factor calculation to be performed by doing a little computation, then a little bit of 
copying and then resuming the computation. This facilitates use of the general approximation approach for matrices 
with a different kind of structure. For example, when considering the problem of incorporating decision feedback into 
the joint detection process, the matrix S has a Toeplitz-block structure instead of btock-Toeplitz structure. It is then 

10 beneficial to periodically alternate between computing and copying. There is no easy way to express this with matrix 
notation alone. Figure 5 provides an illustration of a different form of a Toeplitz-block structure with internal band struc- 
ture and the layout of the elemen ts of R according to whether they should be computed (dark) or copied (light). 
(0037] As a further illustration of the example embodiment of the invention, a particular example of a calculation of 
a matrix factor R is provided in the following paragraphs, in comparison to a calculation using the Cholesky factoring 

is algorithm. 

[0038] The Cholesky algorithm for computing R can be expressed as follows: 



20 



for i in 1 to n 
for j in i to n 



compute re- 



compute ry 



25 




for k in 1 to i -1 



30 



35 



x^x-r h r k j 
if / = i 



else 



=*/r„ 



40 [0039] Note that the division in the case where / * ; can be replaced by a multiplication by pre-computing 1 lr u once 
per row. The interleaving of actual computations and copy operations can then be achieved like this: 
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10 



15 



for i in 1 to n 

for j in i to n 

if we should compute ry 

compute r\j 

else 

copy ry 

copy ry 

if / > * and y > * 

else 



20 



[0040] A numerical example of a correlation matrix S is provided as follows: 



25 



5 = 



4 2 0 
2 4 2 
0 2 4 



30 



[0041) In this example k = 1, that is, the matrix S is not really block-Toeplitz, but rather, it has an ordinary Toeplitz 
structure. The matrix S has a band structure with v = 2. The exact Cholesky factor R, to four decimal places is: 



35 



AO , 



2.0000 1.0000 0.0000 
0.0000 1.7321 1.1547 
0.0000 0.0000 1.6330 



R = 



[0042] As can be seen, R also has band-structure with v = 2 and the elements down the diagonals can be considered 
to be approximately equal. This approximate equality ikjnvien more pronounced for larger matrices. Equation (8) pro- 
vides for the case that only elements r u and r 12 should be computed, and the rest is copied. The approximating algo- 
rithm now starts by computing r u and then r 12 . The element r 13 should be "copied" and is set to zero because the 
element that would be approximately equal to it, r, 12 , does not exist. This is correct because r 13 is known to be zero 
due to the band structure of R. The algorithm then proceeds to copy elements r 22 > r 23 and r 33 from their respective 
approximations r t1 , r 12 and r 22 . Thus, R according to equation (8) is: 




55 



2.0000 1.0000 0.0000 
0.0000 2.0000 1.0000 
0.0000 0.0000 2.0000 



[0043] The calculation of fl< 8 > requires one square regeneration and one division. Using the known approximation 
technique, we would have also had to calculate elen 



eme|it^?^ 
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*< S > = 



2.0000 1.0000 0.0000 
0.0000 1.7321 1.0000 
0.0000 0.0000 1.7321 



requiring two square root operations and two divisions. 
io [0044] Clearty, this is a better approximation than fi< 6 >. But considering that computing the rest of the second row 
does not involve the costly square root operation and division, it would be advantageous to also compute r 23 (but not 
r 33 ). The matrix oriented approach known to those skilled in the art can not express this. The rule for this would be: 

, 5 shall we compute r, t ?s is /SW and) i (V+/-1)fc? (9) 

(0045] This rule computes all non-zero elements of the first kl rows. For our example, with / = 2 the factor matrix 
becomes; 

2.0000 1.0000 o.oooo" 
0.0000 1.7321 1.1547 
0.0000 0.0000 1.7321 



*< 9 > = 



with only two square roots and two divisions. This is a better approximation than tf 5 ), with the same amount of square 
roots and divisions. 

[0046] As will be appreciated by those skilled in the art various modifications may be made to the illustrative embod- 
30 iment without departing from the scope of the present invention. In particular, multi-user detection in a W-COMA system 
can be also effected. In this case data detection from wanted radio signals in the presence of unwanted radio signals 
is expressed according to equation (10): 



35 b dec ^dec(R^y) (10) 

[0047] In equation (10), where R w is the cross-correlation matrix of the system (not a Cholesky factor), and dec () 
is a function which determines a part of the complex plane within which the signals fall. The cross-correlation matrix 
R w , is a Hermitian, positive semi-definite, block-Toepiitz, band matrix. Solving equation (10) efficiently amounts again 
40 to computing the Cholesky factor of As a result of the special structure of R w , as with the correlation matrix A H 
A , the desired Cholesky factor R has the same band structure and Is approximately biock-Toeplitz. 



Claims 

45 

1. A radio communications receiver which operates to detect radio signals and to recover data representative of the 
radio signals in the presence of contemporaneously detected interfering signals, said radio communications re- 
ceiver comprising; 

50 - means to detect said radio signals (6, 8) and to generate digital base band signals representative of said radio 
signals, 

- means to form a correlation matrix (12) from said base band signals, and 

- a data processor (16) which operates to generate a factor (R) of said correlation matrix (S), and to recover 
said data (d) from said factor matrix (R), wherein said data processor (16) operates to generate said factor 

55 matrix (R) by; 

calculating first elements of the factor (R) from the correlation matrix (S) using a factorisation algorithm, in 
accordance with a pre-determined structure of the factor matrix (R) ; and 
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forming seoond elements of the matrix (R) by selectively copying first elements In accordance with said pre- 
determined structure. 



2. A radio communications receiver as claimed in Claim 1, wherein said pre-determined structure of the factor (R), 
5 is determined in accordance with a band structure of the factor matrix (R). 

3. A radio communications receiver as claimed in Claim 1 or 2, wherein the first elements are contained in at least 
one row of said factor (R), and said second elements are contained in at least one second row of said factor (R). 

10 4. A radio communications receiver as claimed in any preceding Claim, wherein said first and second elements are 
present in the same row. 

5./ 
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A radio communications receiver as claimed in Claim 4, wherein the data processor operates to generate said 
factor (R) in accordance with the following steps; 

for each row of said matrix factor (R) (for / in 1 to n); 

• determining, from the predetermined structure; 

• whether the elements of the row (forj 'm i to n) of the matrix (R) should be calculated using said factorisation 
algorithm, and calculating said element using said factorisation algorithm and corresponding elements from 
the correlation matrix (S); or 

- whether the elements of the row (for j in / to n) of the matrix (R) is one of the second elements and can be 
\ copied from said first elements, and copying said elements from said first elements. 

6. A radio communications receiver as claimed in any preceding Claim, wherein the rows are columns of the factor 
25 (R) and the columns are the rows of the matrix (R). 

7. A radio communications receiver as claimed in any preceding Claim, wherein the factor (R) is upper or lower 
triangular. 

30 8. A radio communications receiver as claimed in any preceding Claim, wherein said radio signals and said interfering 
signals are generated in accordance with a code division multiple access scheme. 



9. A radio communications receiver as claimed in any preceding Claim, wherein said factorisation algorithm is the 
Cholesky or Block-Schur algorithm. 

35 

10. A method of detecting a radio signal and recovering data representative of said radio signal in the presence of 
interfering signals, said method comprising the steps of; 

- detecting said radio signals and generating digital base band signals representative of said radio signals, 
*o - [ generating a correlation matrix (S) from said base band signals, 

• V generating a matrix factor (R) of said correlation matrix, and 

• / recovering said data from said factor (R), wherein said factor (R) is generated in accordance with the following 
V operations; 

- /-calculating first elements of the factor (R) from the correlation matrix (S) using a factorisation algorithm, in 
45 I accordance with a predetermined structure of the factor (R) ; and 

- V forming second elements of the matrix (R) by selectively copying elements from said first rows in accordance 

with said pre-determined structure. 



11. A method as claimed in Claim 10, wherein said pre-determined structure of the factor (R), is determined in accord- 
so ance with a band structure of the factor matrix (R). 



12. A method as claimed in Claim 10 or 11, wherein said first elements are contained in at least one first row of the 
factor (R), and said second elements are contained in at least one second row of the matrix factor (R). 



55 13. A method as claimed in Claims 10, 11 or 12, wherein said first and second elements are present in the same row. 



14. A method as claimed in Claim 13, wherein said factor (R) is generated in accordance with the following steps; 
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• for each row of said factor (R) (for i in 1 to n); 
determining, from the predetermined structure; 

- whether the elements of the row (for j in i to n) of the matrix (R) should be calculated using said factorisation 
algorithm, and calculating said element using said factorisation algorithm and elements from the correlation 

5 matrix (S) ; or 

- whether the elements of the row (for j in i to n) of the matrix (R) is one of said second elements and can be 
copied from said first elements, and copying said elements from said first elements. 

15. A method as claimed in any of Claim 10 to 14, wherein the rows are columns of the factor (R) and the columns 
10 are the rows of the matrix (R). 

16. A method as claimed in any of Claims 10 to 15, wherein the factor (R) is upper or lower triangular. 

17. A method as claimed in any of Claims 10 to 16, wherein said radio signals and said interfering signals are generated 
1 3 in accordance with a code division multiple access scheme. 

18. A method as claimed in any of Claims 10 to 17, wherein said factorisation algorithm is a Cholesky algorithm or a 
Block-Schur algorithm. 

20 19. A method of generating at least one factor (R) of a matrix (S) t said factor (R) being upper or lower triangular in 
form, said method comprising the steps of; 

- calculating first elements of the factor (R) from the matrix (S) using a factorisation algorithm, in accordance 
with a pre-determined structure; and 

25 . forming second elements of the matrix (R) by selectively copying elements from said first elements in accord- 

ance with said pre-determined structure. 

20. A method as claimed in Claim 15, wherein said first elements are contained in at least one first row of the matrix 
factor (R), and said second elements are contained in at least one second row of the matrix factor (R). 

30 

21. A method as claimed in Claim 20, wherein said predetermined structure of the factor (R), identifies said at least 
one first row of the factor (R) to be calculated in accordance with the factorisation algorithm, and said second 
elements of the factor matrix (R) in accordance with a band structure. 

35 22. A method as claimed in any of Claims 19 to 21 , and further including the step of; 

- for each row of said factor (R) (for / in 1 to n); 
determining, from the pre-determined structure, 

whether the elements of the row (for j in / to n) of the matrix (R) should be calculated using said factorisation 
*o algorithm, and calculating said element using said factorisation algorithm and elements from the matrix S); or 

- whether the elements of the row (for j in i to n) of the matrix (R) can be copied from said first elements, and 
copying said first elements. 

23. A method as claimed in any of claims 19 to 22, wherein the rows are columns of the factor (R) and the columns 
<5 are the rows of the matrix (R). 

24. A method as claimed in any of claims 1 9 to 23, wherein said radio signals and said interfering signals are generated 
in accordance with a code division multiple access scheme. 

so 25. A data processor which operates to generate at least one factor (R) of a matrix (S), according to the method claimed 
in any of claims 1 9 to 24. 

26. A radio communications receiver as herein before described with reference to the accompanying drawings. 

55 
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